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ABSTRACT 

During the design stage of automatic control systems for 
aircraft, fire control or other physical systems, system com- 
ponents must be selected. Values of the parameters of these 
components must be determined so as to produce the desired 
dynamic properties of the resultant system. 

The method of estimating dynamic characteristics employ- 
ed in this investigation is based on the encirclement criteria of 
Cauchy's Residue Theorem. . A computer was developed to aid in 
this estimation by making certain modifications to the data pre- 
sentation of the Fourier Harmonic Synthesizer. 

This investigation indicated that dynamic characteristics 
may be readily estimated by the method to be described for any 
selected set of system parameters. System parameter specifi- 
cations and limits may also be set when the dynamic character- 
istics are specified. 

Several examples are included to illustrate the techniques 
involved and to indicate the suitability of applying this method to 
aid in the analysis of higher order systems. 

Thesis Supervisor: Sidney Lees 

Title: Assistant Professor 

Aeronautical Engineering 
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OBJECT 



The object of this investigation is to present suitable 
techniques, based on the Cauchy Residue Theorem and utilizing 
a modified Fourier Harmonic Synthesizer as a computer, to aid 
in the estimation of the dynamic characteristics of physical 
systems. 



II 



CHAPTER l 



INTRODUCTION 



In the design of physical systems some of the most important 
considerations are the dynamic properties of the resultant system. 
It will be shown in Summary 1 that the differential equation des- 
cribing the system may be written with the input terms on the 
right side. By setting the right side equal to zero the homogeneous 
equation is obtained. This equation determines thetransient modes 
of the system, which are independent of the input. The actual in- 
put affects the magnitude of the transient modes but has no effect 
on the damping ratios, the oscillation frequencies, or response 
times of the transients. These characteristics are determined 
by the component parameters. Once the system is built, control 
over these dynamic properties is restricted. Since these proper- 
ties have a prominent effect on the performance of the system, the 
engineer must have some means of predicting them during the 
design stage. 

Exact determination of system parameters is often neither 
necessary nor desirable in the early stages of design if limiting 
values of the parameters can be estimated in terms of the dynamic 
characteristics of the system. Such an estimate will enable the 
engineer to set limits on the component specifications in the in- 
itial stages of design. 

The method of estimation of system dynamic characteristics 
described in this thesis is based on a concept first advanced by 
Sidney Lees in Massachusetts Institute of Technology Instrumen- 
tation Laboratory Report R-71 of April 1954 ^ . 

The Mapping Theorem, described in Instrument Engineer- 
ing Volume II ^ ^ and in Theory of Servo Mechanisms, and 
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and 



the Cauchy Residue Theorem described by Guillemin ^ 
others, are the basis for this concept. 

Lees applies this theorem ^ and describes a procedure by 
which dynamic characteristics may be estimated from polynomial 
contours approximated by the computation of certain selected 
points. 

If these contours can be accurately computed and displayed 
by some device, the application of the Cauchy Theorem to this 
problem should be greatly facilitated. 

A Fourier Harmonic Synthesizer built at the Massachusetts 
Institute of Technology Instrumentation Laboratory and described 
by Instrumentation Laboratory Engineering Memo 6445-E-42 
was used for this purpose. Modifications to the existing instal- 
lation and the techniques employed to utilize this device as an 
aid in the estimation of the dynamic characteristics of physical 
systems are described in the succeeding chapters. 
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CHAPTER 2 



MATHEMATICAL BACKGROUND 



The performance of physical systems may often be described 
by linear differential equations with coefficients that are constant 
once they have been decided upon. The properties of the compo- 
nents of the system may be obtained from the equation coefficients. 

The complete solution of the linear differential equation 
consists of the transient solution, plus the forced solution as is 
shown in Summary l. 

The transient solution of the linear differential equation is 
the solution of the homogeneous equation, which is obtained by 
setting the right hand side of the differential equation equal to zero. 



SUMMARY l 



As expressed in Instrument Engineering Volume II ^ , a 
typical ordinary linear differential equation has the form 




n 



m 



(l) 



i=0 



k=0 



The homogeneous equation is 
n 




0 



( 2 ) 



i = 0 



The associated characteristic equation is 
n 




a . z * = 0 



( 3 ) 



i = 0 
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A fundamental theorem of algebra described in Summary 
24-1 Instrument Engineering Vol. II (2) states that every 
equation of nth degree has n and only n roots. The roots of 
Equation (3) are Zy z^, z^, .... z^. In general the roots are 
complex numbers so that 

z i = x i + JYi (4) 

When all roots are distinct, transient solutions of Equation 
(1) have the form 



V 



(tr) 






i=l 



(5) 



Roots of the characteristic equation must have the di- 
mensions of inverse time (from Equation (5) ). Real roots may 
be expressed as 



z. = x. = - 
l l 



(CT). r- 

where (CT)^ = = characteristic time of real root. 

Complex roots may have the form 

z k = x k + + j V 

k k 



( 6 ) 



(7) 



where (DR)^ = damping ratio of kth root, 



u n = undamped natural frequency of kth root, 
k 



16 



and x, = - L , w 
k ^ h ’ 



l 



1 



A 



k vJL/ 

n 

characteristic time of kth root. 



k < CT )k = where (CT)k = Tk 



In polar coordinates, the kth complex root form is 
+ frk = w n, e J 6 = u n, e 3 [* " (ADR) k] 



Z, = X, 

k k 



3 (» - A. ) 
we s k 

n k 



(8) 



where 0 = t. ~ (ADR), = tt - A„ 

k ^k 

and (ADR)^ = A^ = cos~^ = angle of damping ratio of the 
k 

kth root. 



The transient solution determines what oscillatory modes 
can be present, and their mode shapes, and also determines the 
extent to which the dynamic performance of the system will be 
affected T>y the specified input. 

The forced solution reflects the effect of the input and 
determines the relative magnitude of the modes indicated in the 
transient solution superimposed on the forced solution. It also 
determines the system response time due to the effect of the in= 
put on the mode shape but has no effect on the damping or oscil- 
lation frequencies of the second order modes. The response time 
for the system will, in no case, be greater than that indicated by 
the slowest decaying term of the transient solution. 

This investigation is not restricted to any specified input 
and considers only dynamic effects resulting from the roots of the 
characteristic equation. The analysis of this part of the complete 
performance equation gives information on dynamic properties 
which greatly affect the performance of the system for any 
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arbitrary input and it is to this end that the investigation is 
directed. 

The roots of the homogeneous equation expressed as a 
characteristic equation with the variable z are defined in M, I. T. 
Instrumentation Laboratory Report R-71 0) as the dynamic 
characteristics of the associated system. The dynamic charac- 
teristics can be determined if the roots of the characteristic 
equation can be found by factoring or other means. 

Lees (l)discusses the important considerations in the 
investigation of transient solutions of most practical systems. 

One consideration is that the transient solutions decay to zero, 
thus providing absolute stability. In addition, it is necessary 
that the. transient solutions decay rapidly enough to meet specified 
conditions; this is the relative transient stability requirement. 

The damping ratio and the oscillation frequency are also impor- 
tant considerations in investigating the transient solutions. 

These are the dynamic characteristics and are apparent from the 
transient solution or roots of the characteristic equation. In this 
thesis, the terms r, t, and specified by the roots of the charac- 
teristic equation a*re considered to be the dynamic characteristics 
of the physical system. The equivalence of these terms and the 
characteristic equation roots are shown in Figure 1. 




Figure l: Roots of the Characteristic Equation as Represented 

on the Complex Plane. 
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Absolute transient stability of the system exists when all 

of the roots of the characteristic equation have negative real parts, 

x.t x, t 

i. e. , when x^ and x^ are negative, e and e decayjto. zero as 
t — oo . 



Relative transient stability exists when the smallest of the 
magnitudes of the real parts of the roots is greater than a 



specified reference magnitude, i. e. , if 



x. 

l 



mp' 



ref 



then 



(CT). < (CT) r0 . (x^ is most positive real part of all the roots). 

mp 

This has been discussed in greater detail in Instrument Engineering, 
Volume II Chapter 20. 

The damping ratio and oscillation frequency of the system 
are related to the real and imaginary values of the complex root's 
as shown in Equation (7). 

In many practical systems the dynamic characteristics are 
restricted to a certain desired sector or range depending on the 
relative transient stability requirements and specified undamped 
natural frequency and damping ratio limitations. These sectors 
and ranges are shown in Figure 2 below. Specification of a sector 
imposes limitations on the roots and thus the coefficients of the 
characteristic equation, and it is desired to determine the range 
of each coefficient which wilt provide dynamic characteristics 
within specified limits. Figure 2 shows how these limitations may 
be investigated with respect to the characteristic equation complex 
plane root loci. These contours are useful when investigating for 
relative stability, oscillation frequency and damping ratio. In 
most practical systems the specifications impose limitations on 
two or all three of these characteristics and a contour similar to 
Figure 2-d is used. 
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Transient 

Stability 

Limit 



(b) 

Undamped 

Natural 

Frequency 

Limits 



(c) 

Damping 

Ratio 

Limits 




CT 



(d) 

Superimposing Figures 
a, b. and c 

Dynamic Characteristic 
Limits 





Figure 2: Standard Contours 



This investigation is concerned with the techniques con- 
cerning the solution of two problems: determination or estimation 
of the dynamic characteristics, and fixing the limitations and re- 
lationships of the coefficients. Fixation of the coefficients of the 
characteristic equation determines the limits of the component 
parameters in the actual system. 

The simpler case, that of locating the roots of the charac- 
teristic equation with known constant coefficients will be taken up 
first. This case will be treated in some detail since all the 
principles involved are equally applicable in estimating the dy- 
namic characteristics from equations with adjustable coefficients. 

In this presentation, the methods used in estimating the dy- 
namic characteristics of physical systems are based on Cauchy's 
Residue Theorem, as presented by Guillemin (4) and others. 

From Equation (3), 



n 




i = 0 



20 



This is expanded as 



P(z) = a z n + a , z 
v ' n n - l 



2 

+ . . . + a 0 z + a. z + a = 0 



o 



(9) 



From Cauchy's Residue Theorem, if P (z) is a polynomial 
of z, P' (z) is the derivative of the polynomial with respect to 
z, and C is a closed countour on the z complex plane containing 
N roots of P (z) = 0, 



This integral may be interpreted to mean that the argument 
of P(z) increases exactly 2ffN times as the closed contour C is 
traced out in a positive sense, or that the P(z)-plane path encircles 
its origin N times. This is also referred to as the mapping theorepi, 
in Instrument Engineering, Volume II (2), and is illustrated in 
Figure 3. 




C 



( 10 ) 




P(z)T>lane path 
corresponding 
s* to C 



(a) Three roots on z-plane 
within contour 



(b) Three encirclements of the 
origin on the P(z)~plane 



Figure 3. Illustration of Mapping Theorem 
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A closed contour similar to Figure 2-d and corresponding 
to Figure 3-a may be drawn on the characteristic equation root 
complex plane, in which it is desired to investigate for roots of 
the characteristic equation. By substituting values of the 
variable lying along the complex plane contour of Figure 2-d, a 
corresponding contour on the P (z) plane corresponding to Figure 
3-b may be determined. The number of encirclements of the 
F(z) plane origin is the number of polynomial roots within the 
root complex plane contour. As will be shown in the examples, 
the polynomial P (z) of Equation (9) usually includes a constant 
term a Q It is convenient to write the polynomial as P(z)-a Q 
and -a Q becomes the P(z) plane encirclement point. 

The P(z) plane path of Figure 3-b may be calculated by 
hand by substituting values of the variable along the contour of 
Figure 2-d. This becomes impractical since a large number of 
calculations is required if a complete and accurate P(z) plane 
path corresponding to the z plane contour of Figure 2-d is desired. 
Chapter 3 describes how the Fourier Synthesiser may be used to 
obtain the information provided by the P( z) plane path. 
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CHAPTER 3 



APPLICATIONS OF THE FOURIER HARMONIC SYNTHESIZER 



The Fourier Harmonic Synthesizer may be used to facili- 
tate the application of Cauchy's Residue Theorem to this problem, 
reduce the required calculations, increase the accuracy and bring 
other related items of interest to light. 

It has been shown tha 
in the form 

P (z) - a Q = a n R n e 

. . . + a ^ R e 

by setting z = Re ^ . 

This polar form of the equation is suitable for setting into 
the harmonic synthesizer. The method of scaling to obtain full 
scale deflection of the synthesizer and maximum accuracy there- 
by is given and explained in Example l. 

Each of the coefficients (R n , R n ~^) corresponds to the 
radius of the contribution of the particular term of each order as 
the variation is traced out from 0 = 0 to 0 = 2 n. 

The synthesizer traces the polynomial £p (z)Jplane from 
0 = 0 to 0 = 2 7T for the selected R. The plot as recorded on the 
Sanborn Recorder is displayed as separate imaginary (y) and real 
(x) components. The polynomial contour (P(z) plane) which may 
be constructed from the x and y components, is symmetrical 
about the real axis. The real component is a mirror image 
about 180° and the imaginary component crosses through zero 
at 180° with the 180° to 360° portion being an inverted mirror 



: the polynomial in z may be written 



J n0 + a . R "- 1 e 9 + . 

n -l 



j0 



( 11 ) 
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image Therefore, as is indicated in Figure 4, only the 0° to 
180° portion is required. 





(a) Synthesizer Component Traces 

Figure 4 



(b) Polynomial Contour 
Represented by 
Synthesizer Traces 



The P(z) plane may be drawn from the imaginary and real 
components, and the number of encirclements of the constant 
term counted. This, as Cauchy's Residue Theorem states, gives 
the number of roots that may be found at smaller radii (R) than 
that selected. 



The major item of interest is the y (imaginary) cross-over 
(where the y component equals zero) in relation to the constant 
term a Q on the negative real axis. This point is represented by 
a horizontal line at x = - a Q on the x (real) component plot. Only 
the 6 at y cross-overs need be noted from the imaginary com- 
ponent plot, and from the real component plot only the fact of 
whether the corresponding points are more negative or less ne- 
gative than the constant term need be noted in order to sketch 
rapidly and adequately the approximate P(z) plane. If at a y 
(imaginary) cross-over, the x (real) component crosses the con- 
stant term line at that point, roots occur at that radius, one real 
root if 6 = 7r, two complex conjugate roots if 6 equals anything other 
than n. This immediately gives the (ADR) and (DR) of those roots, 
since (ADR) = n-6 and (DR) = cos ( n - 6), as well as the mag- 
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nitude of the roots, since R=u n • The typical P (z) planes of 
Figure 5 will illustrate this. 





Three encirclements of a 



o 



3 roots 



at < l 
n 



(a) P(z) contour for R= l 



One encirclement of a 

o 

1 root at u < . 8 

n 

2 complex roots at u =.8 

r n 

6 = t, = cos j = . 707 



(b) P(z) contour for R = . 8 
Figure 5: Typical P(z) Contours 
The approximate P(z) plane need not be sketched at all 
since the encirclement criteria is immediately apparent from the 
x and y component traces directly. If during the interval between 
two y cross - overs, the real plot crosses the constant line, and 
remains on the opposite side when the y plot again reaches zero, 
one encirclement is indicated. That is, one root on the complex 
plane exists at a radius less than that for which the plot was made. 

This is illustrated in Figure 6. 

0 





One encirclement between 0 and 6 ^ 

No encirclements between 6 ^ and 
©2 or ir 

(a) One Encirclement Indicated by (b) One Encirclement In - 

Synthesizer Traces dicated by P(z) Contour 

Figure 6 
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It should be noted that from the component or p (z) plot of 
just, a single radius, no conclusion should be made as to whether 
the roots indicated are real or complex,, nor estimate the (DR) 
of such roots unless some additional information is known about 
the equation being tested or unless the x plot crosses the constant 
line at or very near the 0 at which the y cross-over occurs. 

Caution should be exercised when attempting to draw concrete 
conclusions in estimating (DR) from plots of x plots crossing the 
constant line near a y cross-over until the effects of R variation 
for the particular equation is fully appreciated and relative mag- 
nitude of synthesizer indications and tape scaling factors am consi- 
dered, Examples (l) and (2) demonstrate, among other things, 
how an error in estimating (DR) may occur by considering only 
the results of a single run 

By making several runs of constant R for varying 0 on the 
machine, a composite plot may be constructed as is shown in 
Example (3). 

The magnitude of the roots may be estimated within the 
limits of the selected R's by the criteria described above. Curves 
connecting points of the various R plots for equal angles (0) pro- 
duce curves that are actually constant 0 curves for varying R. 

The number of encirclements of the constant term by these con- 
stant 0 curves give the number of roots at 0's greater than the 
value of 0 for that curve. This gives the limits on the damping 
ratio, since t, -cos (180°- 0). 

As is indicated in Figure 7, the number of encirclements of 
the constant 0 curves is obtained by considering 0 to be held constant 
until an infinite radius (R) is reached then permitting 0 to complete 
its sweep to 180°, The latter part of the sweep ( @ infinite R) will 
never appear on any trace or plot but must be considered in order 
to determine encirlements. 
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(a) Constant 6 curves for 
Third Order Equation 





(b) Constant 6 curves for 
Fourth Order Equation 



Figure 7: Typical Constant 6 - Variable R curves 
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CHAPTER 4 



TECHNIQUES INVOLVED IN THE TREATMENT OF SPECIFIC 

EXAMPLES 



EXAMPLE 1: Con stant Coefficient Case 

In order to scale the problem in a manner suitable for. 
setting into the synthesizer, the polynomial in z is written in 
the polar form as shown below. 

The polar form coefficients are determined for each of the 
trial radii (R's) selected. The scale setting is then determined 
so that the largest coefficient may correspond to the largest 
synthesizer scale, setting (1. 0) in order that the maximum syn- 
thesizer accuracy may be obtained. This simple set Of slide 
rule ratios lends itself to tabular presentation as is indicated 
in Table I, page 33 which follows the examples. 



P (z) - 52 = z 5 + 13 z 4 + 60 z 3 + 120 z 2 + 124 z 

= R 5 e^ 59 + l3 R 4 e^ 46 + 60R 3 e^ 36 
+ 120 R 2 e j2 6 + 124 Re ^ 

„ j e 

since z = Re 
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EXAMPLE l 



Constant Coefficient Case 



The Sanborn record of the synthesizer component traces are sketched and summarized below for convenience. 
The values at the important points and cross-overs are included and the corresponding P(z) plots sketched in order 
that the effect of R variation may be readily observed. The conclusions and estimates based on the individual runs 
are given with the corresponding traces in order that the technique used in the investigation may be followed con- 
veniently. 



P(z) - 52 = z 5 + 13 z 4 + 60 z 3 + 120 z 2 + 124 z 



y 




5 encirclements; hence 5 roots at R < 6 
At least 1 real root. 

One complex pair may be expected with 
a 6 at 

148° > 0 > 127° 
or 

127° > 9 > 92° 





2 roots @ R = ^2 6 

6 = 145° A ^ * 35° 

(The run for R = 6 indicated the possibility 
of roots @ 6 = 145? The other indication, i.e. 
121^ 0> 92°did not materialize) 

3 encirclements 
3 roots @ R< \f26 

Due to the position of the loop in the P(z) 
plot the only safe estimate is that if there 
is another complex pair, it should occur 
at lOCte 6 < it 



0° 51° 100° 145° tt 





3 roots @ R < 3. 5 

If there is another complex pair at R < 3. 5 
it is more likely to occur between Ul° and it 
then between 59° and 111°. (The basis of this 
estimate is the position of the constant term 
and observing the behavior for many other 
equations) 



0° 59° 111° 
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Example l (cont.) 



0 ° 



110 ° 




l root @ R<1. 2 

(The polynomial contour plot and 
component tracea indicate the 
real root ahould be very close to 
R * l. 2) 





-52 



0 ° 



No roots occur 
at R <. 5 



This checks with previous runs 
since all roots have been located 
at R'a > .5 




The equation 
has roots of 



P(z) = z 5 + I3z 4 + 60z 3 + I20z 2 + I24z + 52 
z= - 5 ± j z *- l± j z = - l 

( u n * 'Jle, A ^ * 11° ) (w n * 'TF, A ^ « 45°) 




The actual values of the roots of the equation agree very well with those 
indicated by the synthesizer. 
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EXAMPLE 2 



Constant Coefficient Case 



Example 2 is of a form very similar to Example l, even to the extent of an increasing magnitude of the 
coefficients with decreasing powers of the variable. 

Example 1: P (z) = z 5 + 13 z 4 + 60 z 3 + 120 z 2 + 124 z + 52 
Example 2: P (z) = z 5 + 5 z 4 + 36 z 3 + 88 z 2 + 108 z + 52 

(The constant term ( a Q * 52 in both examples ) while greatly affecting the roots has no effect on 
the shape or magnitude of the component traces or polynomial contour plot) 

These two examples, while similar in form, show quite different behavior and indicate that an estimate 
of 0 (or A ^ ) based on a single run may be considerably in error. 

P(z) - 52 = z 5 + 5z 4 t 36 z 3 + 88 z 2 + 108 z 




5 encirclements, 5 roots @ R<7 
If complex roots exist, one pair 
may be expected @ 14^ > >0>IIO 
or 11(P>0>84° 







2 roots i> R = V 26 
6-101° A = 79° 

Second estimate from the R = 7 run 
is substantiated. (Compare the 
R = \f26 runs of Example l and 
Example 2) 

3 roots ■§ R < n/26 



3 roots @ R < 3 

If more complex roots exist, it 
cannot be stated positively whether 
they will be at 67$0<136° 
or 136^0 <ir 

The relative magnitude of the trace 
peaka Indicate however, that 0 will 
be closer to 136° than either 67° or n 



2 roots @ R = n/T 
0 = 135° A ; = 45° 

l root (real) @ R < sf2 



0° 52 ° 101 ° 140 ° tr 




32 



Example 2 (cont. ) 



O' 




l root @ R=l 
B-* A^ = 0 

All 5 roots have been located. 

No Indication of roots should be 
found at R’s <1 



ir 




0 ° 





The polonomial: 



P(z) = z 5 + 5z 4 + 36z 3 + 88z 2 +I08z +52 



has roots: 



z= - l± 5 j z*-l±j z = - 1 



^26 , A c = 79°) (u n = -JT , A ; =45°) 



y 




In this example the actual values of the roots agree almost exactly with 
the values indicated by the synthesizer. 
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EXAMPLE 3 



Constant Coefficient Case 



P (z) - 120 = z 3 +l0z 2 + 44 z 
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Example 3 (cont. ) 




2 roots @ R = 4.47 
9 = 116° A^= 64° 




0 ° 





The polynomial 

P(z) = z 3 + 10 z 2 + 44 z + 120 

has root* z = - 6 z = - 2 ± 4j 

( w n ■ 6 ‘- \ - °°) ( w n = 4.47, - 60°) 

The factored form is therefore 

P(s) = (z + 6)(z + 2 + 4j)(z + 2 - 4j) 



The accuracy of B indicated by the synthesizer is very good except when $ approaches 180°. 

Figure 8 shows all the polynomial contours sketched above. It was constructed by taking 
the values from the tapes to provide a series of superimposed polynomial contours for all 
values of R. 

Another series of curves was constructed by connecting points of equal 0 in order that 
constant 6 curves for varying R may be obtained. Cauch's Residue Theorem is applied to 
the constant Q curves. The encirclements are counted and the 6 values determine the sector 
in which roots, if any, are located. 



i 
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Im[P(z)-l20] 
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for Equotion of Third Degree. 



EXAMPLE 3 




Figure 9- Expanded Portion of Region Near Origin in Figure8 for Constant ^-Variable R. 
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Tabulation of Synthesizer Scale Settings for Example 1 



Constant Term 


05 

in 


. 0325 


] . 0621 


1 212 1 


2.27 


1 1 


*'* 1 


00 

tX3* 


Synthesizer Scale Settings 


c/T 


.0442 


o 

03 

r- 

o 


00 

to 


.731 


.862 


o 


1.0 


03 

CO 


. 257 


in 

in 

CO 


.572 | 


1.0 | 


1.0 | 


r- 

to 

03 


00 


CD 

CO 


. 770 


tj* 

O) 

CO 


1.0 | 


.708 | 


.602 


\ W 


03 


CO 


0 T 


o 


.762 


I .217 


.1567 


00 

o 


.0131 


in 

CO 


. 605 


. 393 


.205 


.0236 


. 01445 


to 

o 

00 

o 

o 


. 000505 


Calculation of Coefficients 


124 R 


744 


03 

CO 

to 


434 | 


175.5 | 


149 


*21 


62 


120 R 2 


o 

03 

CO 


o 

03 

CO 


m 

r- 

Tf* 


240 


173 


021 


30 


60 R 3 


o 

CD 

03 

OJ* 


o 

m 

03 

f-* 


o 

00 

m 

03* 


170 


104 


o 

to 


7.5 


4 

13 R* 


O 

CO 

CO 

to* 


o 

03 

r- 

oo* 


m 

to 

03 


52 


VLZ 


CO 


03 

00 


m 

DC 


to 

C- 

o* 


[ OSfr‘8 | 


1 529 


| 5.66 


| 2.50 


- 


, 03125 


Calculation of Powers of R 


DC 


1,296 


to 

c- 

to 


| 151.3 | 


tj* 


2.08 | 


- 


. 0625 


CD 

PC 


216 


to 

03* 

CO 


1 43 


CO 

00 

03* 


1. 73 


- 


.125 


03 

DC 


36 


92 


1 12.3 1 


03 


| !.44 | 




.25 


DC 


to 


in 


in 

CO* 


* 


03 




m 



£ 

u 

o 

H 




to 

dc 



to 

o 

{-. 

03 

3 

O* 

to 



in 

o 






o 

to 



o 

cd 




o 

r- 

t- 



w 

o 



3 

o 



u 

0) 

a 



£ 



w 



CD 



DC 

U 

O 

C*< 



DC 

CO 



o 

to 

05 

of 



n 




o 

CO 
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EXAMPLE 4; Adjustable Coefficient Case 



Example 4 illustrates one analytical method of investigating 
a typical equation with adjustable coefficients in order that the 
limits of the component parameters affecting the coefficient may 
be determined which meet the specified dynamic conditions. This 
method as explained in Example 4 utilizes the algebraic relations 
between the adjustable coefficients and the sums and products of 
the resultant roots. This method is suitable for loweh order 
equations (second or third) with one or two adjustable coefficients. 

EXAMPLE 4 



A typical characteristic equation, of a positional servo- 
mechanism with modified velocity signal damping, is: 



2 ; 



2 7r r 'p 

v _ 

(ps) 1+ 2?r r v 'p 



+ 2 t, 



(PS) ( res 






'p + l =0 



(l) 



where 'p 



w 



n 



(ps) 



Laplace Transform operator 







n 



(ps) 



undamped natural frequency 
of positional servomechansim 



t 



(PS) 



damping ratio of positional 
servomechanism 



r T 
v n 



(PS) 



CT 




non-dimensional character- 
istic time of signal modifier 
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Expanding Equation (l): 



■p 3 + 



2 C 



* (ps) + ^ ^ (ps) (res) + 2 n r 



^ ^ (ps) (res) 2 ^ r v + 'j P + 2 7rr^ 



'P 2 + 



= 0 



( 2 ) 



Or, 



'p^ + (a+c+ b)'p^ + (bc + l)'p + b= 0 



( 3 ) 



where 



a = 2 t 



(Ps) 



b = 



1 



2 7r r 



c = 2 £ 



(ps) (res) 



For zero forced dynamic error, 

FDE = 2 ^ <P S > < res > A , , = 0 for finite u 

“n, , < cm > n (ps) 

( ps ) ' 

where A (angular velocity of controlled member) zjr. 0 



Thus - ^ (ps) (res) 



= c = 0 



and Equation (3) becomes 



'p 2 + (a + b)'p^ + 'p + b = 0 



( 4 ) 



P i 9 

Substituting 'p = — - — > letting p = w e J 

W n, x (r) 

(ps) v 7 
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u n 



and defining R = 



(r) 



U) 



and z = R e 



j e 



n 



(ps) 



3 2 

Equation (4) is: P(z)= z +(a + b)z +z + b = 0 



( 5 ) 



If the roots of the polynomial in z [P(«j| are z ^ , z ^ » and z ^ , 



P (z) = ( z - z : ) ( z - z 2 ) (z - z 3 ) =0 



( 6 ) 



3 2 

z - ( z ^ + z 2 + z ^ ) z + ( z, z^ + z 9 z 9 + z 1 z 9 ) z - 



'l 2 2 3 l 3 



z l z 2 z 3 = 0 



(7) 



For a third order equation with two complex roots, the roots 
will have the form; 



z j = x j - - ct 



Z n = X, 



-2 - ^2 + jy 2 = ' r u n . . + j u n, . 

(r) (O 



V 1 - i „ 2 






z 3 “ x 3 - = - &r u n , , ' l\ , 

(r) (r) 



Using these relations and equating coefficients of like terms in 
Equations (5) and (7) we have: 



a + b = 


a + 2 ^r % , . 

( r ) 


(3) 


1 = 


2 ^r u n, . tt+ w n, . 

(r) (r) 


(9) 


b = 


2 

Of u 

n < i 

(n) 


(10) 
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